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Abstract 

Critical thermoelectric parameters including Seebeck coefficient, electrical conductivity, thermal 
conductivity and figure of merit ZT of one-dimensional coaxial Bi2Te3/Sb2Te3 nanocomposite 
were modeled by following the single carrier pocket and sharp interface assumptions. A calculation 
scheme based on Landauer approach, instead of commonly used Boltzmann transport equation 
(BTE) with relaxation time approximation, was adopted to numerically obtain the transmission 
functions which can be used to evaluate thermoelectric properties. Considerable enhancement of ZT 
was obtained through our modeling and numerical calculation, and the corresponding dependence 
of thermoelectric performance on structure parameters was studied. Finally, the efficiency at the 
maximum power condition for this 1-D system was also investigated. 
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As the demand for renewable energy resources increases dramatically over the last decade, 
the development of high-performance thermoelectric (TE) materials and related devices, 
which can be used to convert waste heat to electricity, has been attracting broader and 
greater interest. 1 Enhancing the thermoelectric figure of merit ZT in bulk materials remains 
difficult due to the inherent inter-couplings between three major physical parameters that 
determine ZT, namely, Seebeck coefficient \S\, electrical conductivity o and thermal con- 
ductivity k. Alternatively, low-dimensional nanostructured materials have been suggested 
to be promising candidates for improving thermoelectric efficiency based upon the proposed 
advantageous properties including the increased phonon scattering due to the high density 
of interfaces for reducing phonon thermal conductivity and the sharp features of electronic 
density of states for improving thermopower, which are absent in their bulk counterparts. 1 
Thus, both theoretical and experimental studies of the thermoelectric properties of low- 
dimensional material systems such as superlattice and quantum dots have been pursued 
extensively. 1-4 

Being one of the efforts to implement 1-D thermoelectric materials with high ZT value 
using low-cost methods, core/shell nanocomposite structures have considerable potentials 
for tailoring electrical and thermal properties for practical device applications. 5-8 However, 
very few theoretical studies of these coaxial-type nanostructures have been reported. In 
the present study, we are herein trying to extend the general framework of thermoelec- 
tric modeling to 1-D core/shell structures to demonstrate optimal thermoelectric figures of 
merit by investigating the microscopic electronic and thermal transport properties of coaxial 
nanocomposites, which can potentially provide better understandings and reliable predic- 
tions of intrinsic properties of thermoelectric nanomaterials. Bi2Te3/Sb2Te3 heterostructure, 
a well-known material system which has already shown outstanding thermoelectric proper- 
ties in its superlattice structure, 9 was chosen for this study. Fig. 1 (a) shows the schematic 
configuration of this coaxial nanocomposite whose total diameter is w, consisting of a Bi2Te3 
core with diameter of d surrounded by a Sb 2 Te3 outer shell. The corresponding simplified 
band structure is also included to qualitatively show the quantum confinement effects that 
result in quantized sub-band levels as carrier pockets. It has been suggested that the inner 
and the outer diameters d and w in SiGe core/shell heterostructures had large impact on 
determining energy separations and the number of sub-band levels that are inside the core 
potential well respectively. 10 For simplicity, therefore, a diameter ratio {d/w) of 1/2 has 
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been used for all the following calculations, which can result in reasonably significant quan- 
tum confinement effects for improving thermoelectric properties. Additionally, single carrier 
pocket and sharp interface assumptions 10 have also been included in electronic transport 
calculations. 

To model thermoelectric properties, TE physical parameters are commonly calculated 
by obtaining the integrals of transport distribution function through solving Boltzmann 
transport equation (BTE) in the relaxation time approximation. 11 In spite of the popularity 
of BTE-related TE modeling strategies, we instead adopted different calculation schemes 
based on Landauer approach which, as suggested by Jeong et al, 12 can provide accurate 
descriptions of TE parameters through simpler implementations. The Landauer expressions 
of basic electronic transport quantities electrical conductivity, Seebeck coefficienct, and the 
electronic thermal conductivity in linear response regime are given by 12-14 

a = 2 -^I j [1/H 
S=^ [V/K] 

where L and A are length and area of a diffusive conductor where electrical and thermal 
conductance G and K e scale with length and can be expressed as electrical and thermal 
conductivity a and k respectively for convenience, and 

with T e (E) = (X e (E) I L)M e (E) defined as the so-called transmission function, 15 X e (E) being 
the mean-free-path for back-scattering of an electron with energy of E, M e (E) being the 
density of electron modes (DOEM) at energy E, 16 and / being the Fermi function in the 
form of f=l/(l+exp((E — £jt)//cbT)). It is worth noting that X e (E)/L is an approximate 
form of X e (E) / (X e (E)+L) in the diffusive limit. Furthermore, DOEM M e (E) in the Landauer 
formalism can be related to the electronic bandstructure and is essentially the product of 
the carrier velocity and the density of states (DOS). For 1-D nanowire case, assuming that 
a single parabolic subband is occupied, the density of electron modes can be evaluated by 17 

M lDe {E) = Q{E-e l ) 
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where G is the unit step function and E\ is the bottom of the first subband. On the other 
hand, the Landauer expression for phonon (lattice) thermal conductivity can be obtained 
through a similar treatment of phonons on a frequency by frequency basis: 18 



Here the photon transmission function can be written as T p (u) = (\ p (u) / L) M p (u) , with X p (u) 
being the mean-free-path for back-scattering of a phonon with a frequency u and M p (oj) 
being the density of phonon modes (DOPM) at a given frequency with a 1-D form of 18 



where ct^m is the lower (upper) frequency limit of the first subband. 

In order to verify the validity of our calculation schemes, a proof-of-concept comparison 
between our calculations and experimental observations are showed in Figure 1(b), where the 
calculated Seebeck coefficients (upper) and phonon thermal conductivities (lower) of n-type 
Bi 2 Te3 and Sb2Te3 bulk materials are plotted along as a function of electrical conductivity 
and temperature respectively with corresponding experimental data from literature. The 
excellent agreements between our modeling and experiments indicate the feasibility of the 
following theoretical modeling. 

Given the methods described above, critical thermoelectric transport coefficients of 
Bi 2 Te3/Sb2Te3 core/sell nanocomposites are modeled at 300 K as a function of geometrical 
dimension - total wire diameter w and relative Fermi level position with respect to the 
conduction (valence) band edge. Figure 2 (a) shows the dependence of figure of merit ZT 
on the size of the coaxial structure and the Fermi energy. It is revealed that a considerably 
high peak ZT value of ~12 can be achieved at a relatively small outer diameter u>~7.5 
nm and with a Fermi level located near valence band edge by p-type doping. However, 
the desired high TE performance with ZT values above 8 was calculated to take place only 
within a very narrow diameter (w) window between 5 nm and 10 nm, and to drop sig- 
nificantly with increasing cross-section size. Further breakdown investigations on different 
component thermoelectric parameters that determine ZT are also showed in Fig. 2(b)-(d) 
to explain the physical mechanism of calculated ZT behaviors. Compared with its bulk 
counterpart, the calculated noticeably higher electrical conductivity o with the peak value 
of ~1.6xl0 7 /f2m and enhanced Seebeck coefficient \S\ up to ~600 uV/K can be obtained 




M lDp {u) = Q(u - Ui)(uf - to) 
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(Fig. 2(b) and (c)), which is not surprising since the proposed sharp features in the electron 
density of states in 1-D systems due to the quantized electron motion perpendicular to the 
potential barrier are expected to produce increased electrical conductivity. The clear size 
dependences of both electrical conductivity and Seebeck coefficient, as shown in the contour 
plots, are believed to arise from the diameter modulated quantum confinement effects on the 
electronic band-structure of studied coaxial structure. Essentially, the separations between 
quantized energy levels and the consequent quantum effects reach their maxima at a small 
diameter value of ~7.5 nm, and then decrease starkly when w further increases, resulting 
in the strong diameter modulation of a and \S\. On the other hand, the Fermi energy 
dependence of shows that a and \S\ values peak near the band edges and decrease when 
the Fermi level Ef further moves into the band. This Ef dependence can be explained by 
realizing the opposite contributions to electrical conductivity and Seebeck coefficient from 
energy levels above and below Fermi level. The density of electron modes on both sides of 
Ef will approach balanced when Ef moves further into the bands, leading to a decreased 
net contribution to thermoelectric transport properties. A similar mechanism also applies 
in size and Fermi level dependence of thermal conductivity, especially the electronic contri- 
bution K e i. The peak value of K e i matches the highest density of electron modes (DOEM), 
as predicted by the expression given above, while the maximum K e i value observed when Ef 
locates in the band-gap can be explained by ambipolar heat conduction mechanism arising 
from the unequal energy flow carried by the balanced electron flow within each band. 

The temperature dependence of the calculated figure of merit and Seebeck coefficient of 
Bi 2 Te 3 /Sb2Te3 coaxial nanostructure has also been studied and plotted respectively along 
with that of bulk materials. The observed strong decreases of Seebeck coefficients (Fig. 3(b)) 
for bulk Bi 2 Te 3 and Sb 2 Te 3 materials when T>300 K are believed to be a result of conduc- 
tivity transition from doped to intrinsic nature at relatively high temperatures, while both 
the excess carriers and quantum confinement shift the peak value of nanostructure Seebeck 
coefficient to higher temperatures with respect to that of bulk materials by suppressing the 
intrinsic conductivity. Moreover, the differences in magnitude of \S\ and peak temperature 
for nanostructures with different wire cross sections and Fermi energy positions apparently 
result from the size and energy level dependence of quantum effect. Similarly, it's not diffi- 
cult to understand the ZT behavior in the presented temperature range, as shown in Fig. 
3(a), by accounting for the temperature dependence of \S\, a and k To predict the practical 
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performance of Bi2Te3/Sb2Te3 core/shell heterostructure, a generic thermoelectric device 
configuration consisting of a cold and a hot electron reservoir obeying a Fermi-Dirac contri- 
bution is considered, as shown in Fig. 4(a), where a symmetrical applied voltage difference 
is assumed for simplicity. 19 The two reservoirs are connected by the coaxial nanostructured 
TE device characterized by its transmission function, and the thermoelectric power will be 
generated when high-energy electrons move from the hot side to the empty states on the 
cold side. Besides, the temperatures for both sides were fixed to be T c =300 K and T H =360 
K (AT/T = 0.2) for efficiency and power calculations in telluride-based ID nanostructures. 
In addition to the thermoelectric efficiency (77), the efficiency at maximum power (rj max p) 
under finite current condition has also been evaluated. The results in Fig. 4(b) show that 
the calculated r/, r] max p and P ma x (scatters) from transmission function approach match the 
solid lines obtained using traditional ZT-based evaluations. 20 All the three quantities in- 
crease with increasing thermoelectric figure of merit, and the structure with the highest ZT 
of ~12 corresponds the highest achievable 77 of 88%*r)c and r\ maxP of 50%*r] C , approach- 
ing the Curzon-Ahlborn limit, the thermodynamically maximum efficiency under maximum 
condition. 

In conclusion, thermoelectric coefficients and figure of merit of coaxial Bi2Te3/Sb2Te3 
nanostructures with circular cross sections have been modeled by studying the corresponding 
transmission functions and related densities of modes for both electrons and phonons based 
on Landauer approach. Strong size, doping and temperature dependences of cr,\S\, k and ZT 
suggest the crucial role of rational material design for reaching the theoretically achievable 
figure of merit up to 12 at 300 K. Additionally, the corresponding thermoelectric efficiency 
and power density of the studied core/shell nanostructure in practical TE devices have also 
been analyzed and r) max p/r)c is found to approach Curzon-Ahlborn efficiency limit of ~51% 
Carnot efficiency when the ZT of Bi 2 Te 3 /Sb2Te3 nanocomposite increases to its peak value. 

The author wishes to thank Dr. Z. H. Qiao and Dr. Bin Wang for help discussions on 
numerical calculations 



6 



Reference 

X L. D. Hicks and M. S. Dresselhaus, Phys. Rev. B 47, 16631 (1993). 
2 Y.-M. Lin, X. Sun, and M.S. Dresselhaus, Phys. Rev. B 62, 4610 (2000). 
3 Y.-M. Lin and M. S. Dresselhaus, Phys. Rev. B 68, 075304 (2003). 

4 A. I. Boukai, Y. Bunimovich, J. Tahir-Kheli, J.-K. Yu, W. A. Goddard III, and J. R. Heath, 
Nature 451, 168 (2008). 

5 K. F. Hsu, S. Loo, F. Guo, W. Chen, J. S. Dyck, C. Uher, T. Hogan, E. K. Polychroniadis, 
and M. G. Kanatzidis, Science 303, 818 (2004) 

6 J. R. Sootsman, R. J. Pcionek, H. Kong, C. Uher, and M. G. Kanatzidis, Chem. Mater. 
18, 4993 (2006). 

7 B. Poudel, Q. Hao, Y. Ma, Y. Lan, A. Minnich, B. Yu, X. Yan, D. Wang, A. Muto, D. 
Vashaee, X. Chen, J. Liu, M. S. Dresselhaus, G. Chen, and Z. F. Ren, Science 320, 634 
(2008). 

8 Y. Ma, Q. Hao, B. Poudel, Y. Lan, B. Yu, D. Wang, G. Chen, and Z. F. Ren, Nano Lett. 
8, 2580 (2008). 

9 R. Venkatasubramanian, E. Siivola, T. Colpitts and B. O'Quinn, Nature 413, 597 (2001). 

10 M. Y. Tang, Master thesis, Massachusetts Institute of Technology, 2004. 

n G. D. Mahan and J. O. Sofo, Proc. Natl Acad. Sci. USA 93, 7436-7439 (1996). 

12 C. Jeong, R. Kim, M. Luisier, S. Datta, and M. Lundstrom, J. Appl. Phys. 107, 023707 

(2010). 

13 U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 
14 C. R. Proetto, Phys. Rev. B 44, 9096 (1991). 

15 S. Datta, Electronic Transport in Mesoscopic System (Cambridge University Press, New 
York, 1997). 

16 S. Datta, Quantum Transport : Atom to Transistor (Cambridge University Press, New 
York, 2005). 

17 R. Kim, S. Datta, and M. S. Lundstrom, J. Appl. Phys. 105, 034506 (2009). 

18 T. Markussen, A. Jauho and M. Brandbyge, Phys. Rev. B 79, 035415 (2009). 

19 N. Nakpathomkun, H. Q. Xu, and H. Linke, Phys. Rev. B 82, 235428 (2010). 

20 G. S. Nolas, J. Sharp, and H. J. Goldsmid, Thermoelectrics: Basic Principles and New 

Materials Developments (Springer, 2001). 



7 




Temperature (K) 

FIG. 1: Schematic diagram of modeled Bi2Te3/Sb2Te3 1-D system: prospective view, cross-section view 
and band diagram, (b) Comparison of Seebeck coefficients as a function of electrical conductivity (upper) 
and phonon thermal conductivity as a function of temperature (lower) for Bi2Te3 and Sb2Te3 respectively 
obtained from experiments and the calculation presented in this study. 
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FIG. 2: Thermoelectric parameters including figure of Merit ZT (a), electrical conductivity a (b), 
Seebeck coefficient \S\ (c) and thermal conductivity k (d) are plotted as functions of the outer diameter of 
coaxial nanostructures and Fermi level energy. Dashed lines represent peak or valley values of 

corresponding parameters 
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FIG. 3: Temperature dependence of figure of merit ZT (a) and Seebeck coefficient |5| (b) for bulk Bi2Te3 
(red solid line), bulk Sb2Te3 (black dashed line), Bi2Te3/Sb2Te3 core/shell with w—20 nm (green solid 
line), w=7.5 nm with £^=-0.05 eV (blue dashed line) and u>=7.5 nm with Ef=-0.l eV (blue solid line). 
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FIG. 4: Schematic diagram set-up of a representative thermoelectric device which can be described by 
its transmission function Tiu(E), as sketched for example, with contact leads acting as a cold and a hot 

electron reservoirs; (b) thermoelectric maximum efficiency (black square), the efficiency at maximum 
power (red circle), both normalized by Carnot efficiency (also labeled in the figure), and maximum power 

density are plotted as functions of figure of merit ZT. 
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